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Abstract 

^ We outline a new method for estimating the cosmic infrared background using the 

spatial and spectral correlation properties of infrared maps. The cosmic infrared background 
f^from galaxies should have a minimum fluctuation of the order of 10% on angular scales probed 
by the DIRBE (Diffuse Infrared Background Experiment) on board COBE (Cosmic Back- 
ground Explorer). We present the theoretical basis for the analysis and discuss the all-sky 
all- wavelengths survey that we conducted. We further show that a linear combination of maps 
O-ftt different wavelengths can greatly reduce the fluctuations produced by foreground stars in 
^ the near-IR, while not eliminating the fluctuations of the background from high redshift galax- 
^ies. The method is potentially very powerful, especially at wavelengths in the mid- to far-IR 
C/3 where the foreground is bright but smooth. 

1. Introduction 



Diffuse background radiation fields contain information about the entire history of the 
universe, including those periods in which no discrete objects exist or can be detected by tele- 
scopic study. The cosmic infrared background (CIB) should contain much of the luminosity 
of the early generations of stars and galaxies, since the cosmic expansion redshifts the wave- 
length of peak luminosity from these sources from the UV and visible bands into the infrared. 
In addition, dust absorbs much of the original UV luminosity and re-emits it in the IR. The 
CIB therefore contains important information about the evolution of the Universe between 
the redshift of last scattering (probed by the microwave background) and today (probed by 
optical surveys). 

Numerous models have been developed to predict the expected properties of the CIB 
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at wavelengths from the near-IR (A ~1-10 //m) to the far-IR (e.g. Bond, Carr and Hogan 
1986; Beichman and Helou 1991; Franceschini et al. 1991). The spectral properties and the 
amplitude of the predicted CIB depend on the various cosmological assumptions like the 
cosmological density parameter, the history of star formation, and the power spectrum of the 
primordial density field among others. A typical prediction over a range of wavelengths probed 
by DIRBE, 1.25-240 fim, is A/a ~ 10 nWm~^sr~^ (e.g. Pagel 1993 and references therein). The 
precise measurement of the CIB levels at various wavelengths would provide information about 
the various cosmological parameters and scenarios. At most wavelengths 10 nWm~^sr~"'^ is 
considerably lower than the residual foreground levels left after the direct foreground modelling 
methods (Hauser 1995a, b), especially in the mid-IR where zodiacal emission is very bright. 
The DIRBE limits from direct modelling (Hauser 1995) have already placed interesting limits 
on the more exotic models of the early Universe. At wavelengths A = (15 — 60)//m there exists 
an interesting possibility of an indirect detection from 7-rays at the level of A/a ^ 6h{^)°^^^ 
nWm~^sr~^ (Dwek and Slavin 1994). Because of the difficulty of direct foreground modelling 
and subtraction, it is also important to try alternative methods as well. 

One such method can be based on the angular correlation function analysis of the 
DIRBE maps, the groundwork for which has been laid in Kashlinsky, Mather, Odenwald and 
Hauser (1995, 1996; hereafter Paper I, Paper II). Here we review the theoretical background 
for the method and report on the limits it sets on the CIB. 

2. Theory 

The correlation signal for the CIB has been studied for a variety of cosmological models 
(e.g. Bond, Carr and Hogan 1986, Wang 1991, Coles, Treyer and Silk 1990, Paper I). In general 
the intrinsic correlation function of the diffuse background, C'(6) = {X6I\(x) ■ X6I\(x + 6)) 
with 61 = /—(/), produced by a population of emitters (e.g. galaxies) clustered with a 
3-dimensional correlation function ^(r) is given in the small angle limit (6 << 1) by: 

Cie) = Ae{z){—^f[^\z){l + zfVTT^,]dz (1) 

where it was assumed that the cosmological constant is zero, ^{z) is the factor accounting for 
the evolution of the clustering pattern, and 



Ae{z) = 2R-^' a\Jv' + 0^ )dv (2) 

Here Rh = cHq^ , and x(z) is the comoving distance. After convolving (1) with the beam the 
zero-lag signal for the DIRBE sky becomes 

Cm= r AAi)i'-^fl*Hi)il + :fVT+n^]<li. (3) 

For the top-hat beam of DIRBE the window function is W(x) = [2Ji(x)/x]^ and «?=0.46°. In 
the above equation Ps^ik) is the spectrum of galaxy clustering at the present epoch. If ^ is 
known the measurement of C(6') can give information on the diffuse background due to material 
clustered like galaxies (e.g. Gunn 1965, Peebles 1980). Such a method was successfully applied 
in the V (Shectman 1973,1974) and UV bands (Martin and Bowyer 1989). In the rest of the 
paper we omit the subscript i9 in (3) with C(0) referring to the DIRBE convolved zero-lag 
correlation signal, or the variance of the map {(\6I\)^). 
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As was shown in Paper I for scales dominating the integral in (2) the 2-point correlation 
function at the present epoch can be approximated as a power-law. In that case in the 
Friedman- Robertson- Walker Universe A^)(z) would reach a minimum at z > 1 whose value 
is almost independent of SI (Paper I). This argument was generalized in Paper II to a more 
realistic ^(r) which was computed using the correlation data from the APM (Maddox et al. 
1990) survey with the APM power spectrum taken from 1) the inversion technique of Baugh 
and Efstathiou (1993) and 2) the empirical fit to the APM data on the projected galaxy 
correlation function from Kashlinsky (1992). Thus we can rewrite (3) as an inequality in order 
to derive an upper limit on a measure of the CIB fiux from clustered material from any upper 
limit on C(0) derived from the DIRBE data: 

(A/a).,™. < S^CXO) (5) 

where B = 1/ y/mm{A,}{z)} = (11 - 13) over the entire range of parameters, and the measure 
of the CIB fiux used in (5) is defined as [iXIx)z,rm.s? = J{^f[^^iz)i'^ + zf^J\^^z\dz. 
The latter is ~ \{^^fdz since the term in the brackets has little variation with z for two 
extremes of clustering evolution, when it is stable in either proper or comoving coordinates 
(Peebles 1980). Thus the CIB produced by objects clustered like galaxies should have sig- 
nificant fiuctuations, ~ 10% of the total fiux, on the angular scale subtended by the DIRBE 
beam. In general, as eq. (3) shows if the bulk of the CIB comes from higher redshift this 
would lead to smaller relative bljl and vice versa (cf. Wang 1991). 

3. Data and analysis 

The COBE was launched in 1989, and included the DIRBE instrument with the aim of 
studying the IR emission in ten wavebands from the near-IR (J,K,L bands at A=1.25, 2.2 and 
3.5 pim) to the far-IR extending to A=240 //m. Four of the bands were designed to match the 
IRAS. The instrument is a photometer with an instantaneous field of view of 0.7x0.7 degrees. 
It was designed to have instrumental noise less than 1 nWm~^sr~^ in the first eight bands; 
at the two longest wavelengths the instrumental noise is considerably higher and this would 
make, as discussed below, the correlation method less efficient there. We have used the map 
derived from the entire 41 week DIRBE data set available from the NSSDC. A parameterized 
model developed by the DIRBE team (Reach et &\. 1995) was used to remove the time- varying 
zodiacal component from each weekly map (see Paper I for more details). The maps were 
constructed for each of the ten DIRBE wavebands. The DIRBE maps were pixelized using 
the quadrilateralized spherical cube as described by Chan and O'Neill (1974) and Chan and 
Laubscher (1976). 

After the maps were constructed the sky was divided into 384 patches of 32x32 pixels 
or ~ 10°xlO° each. Each field was cleaned of bright sources by the program developed by 
the DIRBE team. The fiux distribution in each patch was modelled with a smooth fourth 
order polynomial component. Pixels with fiuxes > #cut standard deviations above the fitted 
model were removed along with the surrounding 8 pixels. Three values of #cut were used: 
= 7,5,3.5. Since any large-scale gradients in the emission are clearly due to the local 
foregrounds a third order polynomial was removed from each patch after desourcing. 

Three types of foregrounds provide dominant contributions at various bands: in the 
near-IR bands (1.25-4.9 //m) the main foreground contribution to the zero-lag signal comes 
from the Galactic stars, in the mid-IR (bands 12-60 //m) from the zodiacal light emission, and 
in the far-IR bands (100-240 //m) the main contribution is from the cirrus clouds in the Galaxy. 
The structure of the near-IR maps and the distribution of fiuctuations within the various fields 
have been discussed at length in Paper I; the histograms for these bands are shown there in 
Fig. 5. They are highly asymmetric with respect to blx and are consistent with those from 
simulated stellar catalogs. The distribution at mid- to far-IR bands is more symmetric and 
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comes from a more extended emission. Furthermore, because of the extended character of the 
foreground emission at mid- to far-IR the decrease of the width of the histogram (another 
measure of C(0)) with #cut is significantly less prominent than for the near-IR where stars 
provide the dominant foreground. Removal of the large-scale gradients does not lead to any 
noticeable difference in the histograms at the near-IR bands. On the other hand, at longer 
wavelengths the bulk of the foreground comes from extended emission, and removing large- 
scale gradients, which are clearly of local origin, changes the structure significantly. For the 
same reason there is more noticeable decrease in the measured C(0) with #cut after removing 
the local gradients. 

In the near-IR where Galactic stars dominate the fiuctuations of the foreground there 
is significant decrease in the minimal value of C(0) with #cut- The reason for this is that the 
structure of the foreground is dominated by fiuctuations in point sources (stars) affecting the 
small scale structure of the correlation function. In the near-IR the histogram distribution of 
C(0) on the sky becomes increasingly asymmetric with decreasing #cut and the bins with low 
C(0) patches at higher clipping thresholds now have essentially the same C(0). On the other 
hand, at mid- to far-IR there is very little change in the distribution of the zero-lag signal with 
decreasing #cut ; this is due to the fact that the foreground emission is extended with little 
small-scale structure. The latter in turn leads to much smaller upper limits on the value of 
C(0) at mid- to far-IR. At bands 9,10 the instrumental noise is substantially larger because 
bolometers are used, and the correlations analysis of the DIRBE maps leads to significantly 
less interesting limits at these wavelengths (140, 240 //m). 

The C(0) signal is very anisotropic, showing the shapes of the Galactic and Ecliptic 
planes, and hence most of it comes from the foregrounds. In fact, in the near-IR bands where 
Galactic stars dominate the foreground there is a strong correlation between the value of 
C(0) and the Galactic latitude. At longer wavelengths the correlation with h is significantly 
less pronounced due to appreciable (and at times dominant) contribution from the zodiacal 
foreground. The map of C(0) on the sky clearly shows that the plane of the Galaxy dominates 
the maps at the near-IR (1.25-4.9 /um), although at 4.9 //m the ecliptic plane would dominate 
the fiux maps. At the mid-IR bands, where zodiacal emission is significant, the C(0) map 
shows both the Galactic and the ecliptic plane. At the far-IR where most of the emission 
comes from the cirrus, the Galactic plane again becomes prominent. 

4. Colour subtraction method for reducing foreground 

The signal from the foreground radiation may also be disentangled from that of galax- 
ies by considering the properties of the foreground stellar emission that are different from 
the GIB. The extragalactic objects have somewhat different colors than the Galactic stars, 
and they are redshifted as well. Furthermore, galaxies of different morphological types have 
different stellar populations and this too should help distinguish the colors of the foreground 
emission from our Galaxy (type Sc) from those of more representative galaxies (E/SO, Sa, etc.) 
which are expected to dominate the GIB at J, K and _L-bands. 

We applied these ideas so far only to the three near-IR DIRBE bands. The near-IR 
DIRBE maps exhibit a clear color with a relatively small dispersion. We define the color 
between any two bands, 1 and 2, as ai2 = h/h- Within each field the dispersion of colors 
is small, between 5% and 10%. The colors confirm that most of the light of the Galaxy in 
these bands comes from K and M giants (Arendt et al. 1994). This empirical property of the 
constancy of color found in the DIRBE near-IR maps can be used to reduce the foreground 
contribution to C(0) by taking the scaled difference of two maps. On the other hand, as 
discussed above most GIB in the near-IR is expected to come from galaxies whose spectra are 
significantly redshifted and thus the GIB is not expected to be washed out by this procedure. 
These observations are the basis for the color subtraction method we present below. 

There are two independent contributions to the measured fiuctuations in the linear 
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combinations of the DIRBE maps: the one that comes from the <o<a/ foreground (stars, zodiacal 
light, and cirrus), which we denote as Ch,_a(0), and the one that comes from the extragalactic 
background which we label Cj,a(0). We use the A notation to indicate that these variances 
refer to the differences of maps at two different wavelengths. The total Ca(0) is the sum of 
these two components. Let us assume that the Galactic and zodiacal foreground component 
has the mean color between any two wavelengths 1 and 2 of = (a) and the dispersion in 
the color is = {{^'^)^)y where 6a is the deviation of the pixel color from the mean of the 
pixels in the field. Assuming that the color fiuctuations are uncorrelated with the intensities, 
the zero-lag signal due to the foreground in Band 1, Ch,_i(0), can be written as 



Note that the approximations used in this equation are not required for the validity of the 
analysis of the maps, but are only to indicate that major improvements are achievable with 
color subtraction. Let us construct the quantity which depends on the parameters in two 
bands (we define 6i = li — (li)), Ah,_i2 = — P^*,2, where the quantity [3 is for now an 
arbitrary number which will be quantified later. The zero-lag signal in the "color-subtracted" 
quantity defined above would be given by 



It is clear from the above that if we were to choose [3 close to the mean color for the foreground, 
[3 ~ 512, then the largest term in (7) could be made nearly zero and the contribution of the 
Galactic foreground confusion noise to the zero-lag signal Ch,_a(0) in the difference map would 
be minimized. Furthermore, ii s <^ S, that contribution can be decreased over the one specified 
in (6) by a fairly substantial amount. 

Measurement shows that indeed we have s <^ S. Thus the variance in the color- 
subtracted fiuctuation A due to the confusion noise can be reduced by a substantial factor, 
with a minimum when [3 = ^12. The plots of Ca(0) for the data from all five fields as a function 
of [3 for J — K, J — L and K — L subtraction, for maps of various regions selected from sky both 
before and after zodiacal light removal show that the minimum of Ca(0) is indeed quite deep 
(see Paper I for more details). Therefore, in constructing the color-subtracted quantities we 
use the values of /3 for each field that correspond to the minimum of Ca(0). Unlike the single 
band zero-lag signal, which would correspond to /3 = 0, the minimum of Ca(0) is distributed 
significantly less anisotropically. This in turn leads to interesting and low limits on the near-IR 
GIB as discussed in Paper L 

5. Conclusions 

As a result of this survey we obtained interesting limits on the GIB from the C(0) mea- 
surements. These results with the relevant discussion will be presented in Paper II. Assuming 
the "canonical" number of 10% fiuctuations on the scale of the DIRBE beam the results lead 
to interesting upper limits on the GIB fiuxes from clustered material, particularly at mid-IR 
bands where the foreground is bright but smooth. Because of DIRBE detector noise, they 
are less interesting at 140 and 240 //m. At the J,K,L bands the no-evolution (Paper I) or 
reasonable evolution (Veeraraghavan 1995) models of the zero lag signal in the GIB produced 
by normal galaxies give numbers around (1.5-3, 0.5-1.5, 0.5-1.5) nWm~^sr~^ respectively. Our 
limits in band 1 (1.25 //m) are almost an order of magnitude higher, but in bands 2, 3 (2.2, 
3.5 //m) they are comparable in magnitude. If there was an extra population of galaxies at 
early times they would contribute an additional signal (cf. e.g. Gole, Treyer and Silk 1992). 
Also any other sources of emission from clustered material in the early Universe in addition 
to normal galaxies would increase the above theoretical estimates. 



a, 1(0) ~ {Sl^ + s?2)C^*,2(0) + sl^ih^^y. 



(6) 



a,A(0) = (A2_i2) = [(5i2 - f^f + S?2]C?*,2(0) + S?2(^*,2)^ 



(7) 
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The color subtraction works well to suppress the fluctuations due to the Galactic 
foreground. On the other hand the predicted extragalactic background fluctuations are unlikely 
to be removed by the color subtraction method because the redshifts give the background a 
difli"erent color than the foreground. The extragalactic signal would change by a different 
amount, and may even increase because of the combination of redshift and morphological 
effects. This is discussed in more detail in our Paper I. Since the Galaxy colors are typical 
of general stellar populations our method effectively reduces the contribution from the nearby 
galaxies, but because of the redshift effects on galaxies the overall zero-lag signal remains 
at the levels of the single band cases. Thus the color subtracted signal probes galaxies in a 
window of redshifts centered on an early epoch; the distant galaxies do not have a single color, 
and can not be removed in the process. 
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